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NORMAL SHIFT IN GENERAL LAGRANGIAN DYNAMICS. 


Ruslan A. Sharipov 


Abstract. It is well known that Lagrangian dynamical systems naturally arise in 
describing wave front dynamics in the limit of short waves (which is called pseudoclas- 
sical limit or limit of geometrical optics). Wave fronts are the surfaces of constant 
phase, their points move along lines which are called rays. In non-homogeneous 
anisotropic media rays are not straight lines. Their shape is determined by modified 
Lagrange equations. An important observation is that for most usual cases propa¬ 
gating wave fronts are perpendicular to rays in the sense of some Riemannian metric. 
This happens when Lagrange function is quadratic with respect to components of 
velocity vector. The goal of paper is to study how this property transforms for the 
case of general (non-quadratic) Lagrange function. 


1. A SIMPLE PRELIMINARY EXAMPLE. 

Description of most wave phenomena is based on wave equation. This is second 
order partial differential equation of the following form: 


_ V 

0,2 2^ 


dt 


ay 

9a;* ^ 


= 0 . 


( 1 . 1 ) 


Here t is time variable, while x^, , and x^ are spatial Cartesian coordinates. 

Parameter c in first term is the velocity of wave process described by the equation 

(1.1) . This is sound velocity for sound waves in gases, liquids, or solid materials, and 
this is light velocity for light waves in refracting media. For homogeneous media c is 
constant, but below we consider non-homogeneous media, where c = c{t, x^, x^, x^). 

Function ijj = exp(i (to t — ki x^ — k 2 x'^ — k^ x^)) is a solution of wave equation 

(1.1) for the case c = const. It describes a plane wave. Here w is a frequency of 
wave, while fci, fc 2 , and fca are components of wave vector k. Frequency oj and wave 
vector k are related with each other as follows: 


w = c ■ |k|. (1-2) 

The relationship (1.2) is called dispersion law. Short wave limit corresponds to 
the case of high frequency, when a; —> oo. Below we consider this case for non- 
homogeneous media with c const. Therefore we cannot use simple exponential 
solution tp = exp(i {uj t — ki x^ — /c 2 x^ — k^ x^)) of wave equation (1.1). However, 
we can look for the exponential solution with large parameter A —> oo: 


tP = 


OO 

E 

a=0 


(zA)“ 


j\s 


(1.3) 
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Substituting (1.3) into the equation (1.1), we get the following equation for S: 



(1.4) 


This is well-known eikonal equation (see Chapter VII in [1]). Suppose that refracting 
properties of medium do not change in time. Then c = c{x^,x'^,x^). In this case 
we can consider a wave with constant frequency lo = X. For such wave function S 
in eikonal equation (1.4) is taken to be linear function in time variable t: 

S = t — ip{x^ ,x'^,x^). (1-5) 


Eikonal equation (1.4) then is written as the equation for gradient of (p: 


E 


dtp 

dx^ 


-- = 0- 


( 1 . 6 ) 


Let’s denote —l/(2c^) = U and then let’s write the equation (1.6) as 

3 


E -L {7(x^a:^a;^) = 0. 


(1.7) 


Here S/ip, X/ 2 P, and V 3 P are components of gradient Vp. We denote it by p and 
treat as a vector field in tree-dimensional space 


P = V(/9 = 


dpjdx^ 

dpjdx^ 

dpjdx^ 


( 1 . 8 ) 


If we substitute components of vector (1.8) into (1.7), then we can write (1.7) as 

i7(Vi(,?, V2i^, V3(^,a:\a;^,a;^) = 0, (1.9) 

where function H = H{pi,p 2 ,P 3 ,x^,x^,x^) looks like Hamilton function of a par¬ 
ticle of unit mass m = 1 in potential field U = U{x^, x"^, x^): 


^ = E 


M! 

2 


+ u. 


( 1 . 10 ) 


Let p = p{x^ .x^ ,x‘^) be a solution of the equation (1.7) and let p = Vp be 
corresponding momentum vector field (1.8). Let’s consider integral curves of vector 
field p. They form two-parametric family of curves in 


{ x^ = 

x^ = x^(t,y\y^), 
x^ = x^(t,y\y‘^) 


( 1 . 11 ) 
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defined by solutions of the following system of ordinary differential equations: 


X^=Pl, i'^=P2, x^=P3- (1-12) 

Differential equations (1-12) can be written as Hamilton equations: 

f) M 

i:* = —, i=l, 2, 3. (1.13) 

UPi 

Now let’s calculate time derivative p for the momentum vector (1.8) due to the 
dynamics determined by differential equations (1.12): 

. ^ ^ -k -k ^ 

* ^ dx^ ^ dx'- dx^ ^ dpk dx^ dx'^ ’ 

k—1 k—1 k—1 ^ 

Remember that (/? is a solution of the equation (1.9). Differentiating (1.9), we get 


dH{Viip,V2P:^3P,x'^,x"^,x^) 

dx^ 


^ dH 
^ 9Pk 

k—1 


dx^dx^ 



Comparing the above two equalities, we derive differential equations 


Pi = - 


dH 
dx^ ’ 


z = 1, 2, 3. 


Both (1.13) and (1.14) form complete system of Hamilton equations 


(1.14) 


X 


I 


dH 
dpr ’ 



(1.15) 


with Hamilton function (1.10). 

Note that Hamilton equations (1.15) is a system of 6 first order ODE’s. Its 
solutions define five-parametric family of curves in K.^: 


( x^ = x\t,y^,y'^,y^,y*,y^), 

< x"^ = x'^{t,y^,y‘^,y^,y*,y^), (1.16) 

[ = x^(t,y^,y‘^,y^,y^,y^). 


Curves (1.11) form two-parametric subfamily in five-parametric family of curves 
(1.16). They are distinguished by the following two properties: 

1) curves (1.11) correspond to zero level of energy H = 0; 

2) curves (1.11) are perpendicular to level surfaces of the function (p{x^,x‘^, x^). 
First property follows from (1.9). Second is obvious, since curves (1.11) are directed 
along gradient vector (1.8). One can calculate complete derivative of the function 
(p{x^, x'^, x^) with respect to parameter t along these curves: 


dp 

dt 


n = Yl 



( 1 . 17 ) 
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Curves (1-16) defined by Hamilton equations (1-15) and restricted by the above 
conditions 1) and 2) are called characteristic lines for nonlinear first order partial 
differential equation (1.9). They are used in order to construct solutions of this 
equation as described just below (see also [2] and [3]). 

Let’s take some smooth surface ct in We assume that cr is level surface 
with (^ = 0 for the solution ip{x^, , x^) of the equation (1.9) that we are going to 

construct. Denote by and inner curvilinear coordinates of points on cr. Then 
we can write the equations determining points of a in parametric form: 

( x^ = x\y\y‘^), 

\x^=x\y\y% (1.18) 

[ x^ = x^{y'^,y'^). 

At each point of a we have unit normal vector n. Let’s denote it by n = n{y^,y'^). 
Assuming a to be orientable, we can take n(y^, y^) to be smooth function of y^ and 
y^. Under these assumptions we define vector function 

p = ^ • n (1-19) 

on cr, getting scalar factor v = v{y^,y'^) from the following equality: 

H{vni, V 712, V n^,x^ ,x^ ,x^) = 0. (1-20) 

Then we use vector function (1.19) in order to set up Cauchy problem 

=x\y^,y^), p^\ =p,{y^,y^) (1.21) 

lt=0 li=0 

for Hamilton equations (1.15). Solving this Cauchy problem (1.21), we obtain 
two-parametric family of characteristic lines given by functions (l.H) that extend 
initial functions (1.18). They possess property 1), since we determine v by (1.20). 
They also possess property 2), since we determine p by (1.19) (at least for initial 
surface cr). These characteristic lines fill some neighborhood of initial surface cr. 
Therefore we can treat t, y^, y^ as curvilinear coordinates in and consider (1.11) 
as transition functions to these curvilinear coordinates. Then integral 

t 

ip = J ndt, (1.22) 

0 

where H is given by right hand side of (1.17), yields a solution of partial differential 
equation (1.9) expressed in curvilinear coordinates t, y^, and y^. This solution 
(1.22) satisfies zero boundary-value condition on cr: 

ip\ =0. (1.23) 

1(7 

In other words, (1.23) means that a is zero level surface for the function (p. 

Note that in curvilinear coordinates t, y^, y^ initial surface a is given by the 
equation t = 0. However, other level surfaces of the function p are not given by the 
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equations t = const. In order to change this situation we should choose another 
set of curvilinear coordinates s, y^, where s = Lp{t,y^,y‘^). This means, that 
we change parametrization of characteristic lines (1-11) without changing them as 
geometric sets of points. In new parameter s characteristic lines of the equation 
(1.9) are given by modified Hamilton equations 


1 dH 
n dpi ’ 


, _ 1 dH 

ft dx'‘ 


(1.24) 


where denominator fl is determined by right hand side of (1.17). In our particular 
case, when function H is given by formula (1.10), we have O = (pi)^ + (p 2 )^ + (ps)^- 
Hence 17 0 for p 0. 

Let’s fix new curvilinear coordinates s, y^, y^. Here (f{s, y^, y^) = s by definition. 
Now let’s return to initial wave equation (1.1) and to formula (1.5) for the function 
S in asymptotical power expansion (1.3). It’s important to note that t in (1.5) do 
not coincide with t in (1.11) and in Hamilton equations (1.15), where t was used 
as a parameter on characteristic lines of the equation (1.6). Therefore now in the 
expression for S we have both s and t (and t is time variable again): 


S = S{t,s,y'^,y‘^) = t- s. 


For exponential factor in (1.3), taking into account that A = w, we get: 


_ giaj(i-s) 


(1.25) 


Right hand side of (1.25) corresponds to plane wave propagating in the direction of 
s-axis. In original Cartesian coordinates x^, a;^, x^ this looks like non-plain wave 
propagating along characteristic lines of the equation (1.9). Level surfaces of the 
function p are the surfaces of constant phase in such wave. They are called wave 
fronts. The equation t — s = const, when transformed to Cartesian coordinates 
x^, x^, x^, describes moving surface, that gradually passes positions of level surfaces 
of the function p. This process is called wave front dynamics. It’s very important 
that this process can be understood as a motion of separate points of wave 
front, each obeying modified Hamilton equations (1.24). For this reason these 
equations are called the equations of wave front dynamics. 

Another important point concerning wave front dynamics, that we noted above, 
is that level surfaces of the function p are perpendicular to characteristic lines 
(1.11). Therefore wave front dynamics is a normal displacement (or normal shift) 
of initial surface cr along trajectories of modified Hamiltonian dynamical system. 


2. More complicated example. 

Let M be some Riemannian manifold. Denote by V standard covariant differen¬ 
tiation determined by metric connection F in M. The following differential operator 
is called Laplace-Beltrami operator in the manifold M : 


*=i i=i 


( 2 . 1 ) 
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Here are components of metric tensor in local coordinates x^, ... , a;", while 
Vi and Vj are symbols of covariant derivatives in these local coordinates. Differ¬ 
ential operator H is called fiberwise spherically symmetric if it is represented as a 
polynomial of Laplace-Beltrami operator (2.1): 

m 

H{p,D) = Y,au{p)AK ( 2 . 2 ) 

fc =0 

Here p is a point of M and D is a formal symbol for differentiation. Coefficients 
ao, ... , Qra in (2.2) are arbitrary smooth functions of p G M. Note that H is scalar 
operator. Operator (2.2) can be applied either to scalar field or tensorial field in 
M, yielding the field of the same type as that it was applied to. 

Now let’s add differentiation in time variable dt = d/dt and let’s introduce large 
parameter A to (2.2). As a result we get differential operator 


H{p,X-^D)=J2J2 


s—0 k—0 


ask (p) 
{i A)®+2'= 


a® A'' 


The following differential equation in M is an analog of wave equation (1.1): 


H{p, = 0. 


(2.3) 


Short wave asymptotics X ^ oo for this equations is described by the same asymp¬ 
totical expansion (1.3) as in case of standard wave equation. Coefficients ask(p) in 
(2.3) do not depend on t. Therefore we can choose S to be linear function of t: 

S = t- p(p), (2.4) 


just like it was in (1.5). Substituting (2.4) into (1.3) and substituting (1.3) into 
(2.3), we derive differential equation for phase function (p(p) in (2.4): 

m / m \ m 

E E = E = 0- (2.5) 

k—0 \ s—0 / k—0 

Here |V(p| is modulus of covector field Vp measured in Riemannian metric g: 

n n 

j=i 


Let’s denote Vp by p as it was done above in section 1 (see formula (1.8)): 

difldx^ 

p = Vp = : 

dipjdx^ 

Now we can write (2.5) as polynomial equation with respect to components of p: 


E^'=(p) ■ IpI^^ = 


(2.6) 
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Let’s denote left hand side of (2.6) hy H = H{p,p). The equation (2.5), written as 

Hip,Vip) = 0, (2.7) 

is exact analog of the equation (1.9) from section 1. Further steps in solving this 
equation are quite similar to those in section 1 (they are described in details in 
paper [4]). Below we shall not discuss them. However, we shall point out most 
important features of wave front dynamics in the limit of short waves A = w —> oo 
for generalized wave equation (2.3). They are the following ones: 

— wave fronts are level hypersurfaces at = {p € M : pip) = t} for the function 
p{p), where p{p) is a solution of differential equation (2.7); 

— time evolution of wave fronts at in M can be described in terms of motion of 
their points obeying modihed Hamilton equations 


1 dH 
H dpi ’ 


. _ 1 dH 

H dx'‘ 


where Hamilton function H is determined by left hand side of (2.6) 


H{p,p) = ^bk{p) ■ IpP'”, 

and denominator Q in (2.8) is determined by formula 





( 2 . 8 ) 


(2.9) 


— wave front dynamics for wave equation (2.3) in short wave limit A ^ oo is a 
normal shift of of initial wave front hypersurface cr along trajectories of modihed 
Hamiltonian dynamical system (2.8), this means that orthogonality of wave 
fronts at and trajectories of shift is preserved in time; 

— normal shift of hypersurface cr is initiated by Cauchy problem data 


xM =x"{p), pi\ =i^{p)-ni{p) (2.10) 

li=0 lt=0 

for the equations (2.8), where ni{p) are covariant components of normal vector 
n(p) at the point p G a, and ^{p) is a scalar factor determined by the equation 

H{p,v-nip)) =0. (2.11) 


For us the most important feature of wave front dynamics, among those listed 
above, is the phenomenon of normal shift. It was revealed in simplest case con¬ 
sidered in section 1. It is also present in more complicated case related to some 
Riemannian metric. Our aim below is to reveal this phenomenon for the case of 
general Hamilton function H, which is not restricted by formula (2.9). In order to 
do this we need to introduce geometrical technique, which is not new, but never¬ 
theless, is not commonly known. It seems to me, that this technique first appeared 
in Finslerian geometry (see [5] and [6]). We used this technique in [7-22], where 
theory of Newtonian dynamical systems admitting normal shift was developed (see 
also theses [23], [24], and recent papers [25-29]). 
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3. Extended tensor fields. 

Let’s consider Hamilton function (2.9). It depends on two arguments p and p, 
where p is a point of manifold M, while p is cotangent vector at the point p, i. e. 
p is an element of cotangent space T*{M). Both p and p, taken together, form a 
pair q = (p, p) which is a point of cotangent bundle T*M. This means that id is a 
scalar field in cotangent bundle T*M. But we shall treat it as extended scalar field 
in M as defined below. Let’s consider the following tensor product: 


r times 

t;(p, m) = 'rp(M)® ® t;(m) 0... ® t;(m) 

'-V-" 

s times 


Tensor product TJ(p, M) is known as a space of (r, s)-tensors at the point p G M. 
Pair of integer numbers (r, s) determines type of tensors. Elements of TJ(p, M) are 
called r-times contravariant and s-times covariant tensors or simply (r, s)-tensors. 

Definition 3.1. Extended tensor field X of type (r, s) in M is a tensor-valued 
function that maps each point q = (p, p) of some domain G C T*M to a tensor of 
the space TJ(p, M). If G = T*M, then X is called global extended tensor field. 

Note a trick: arguments of extended tensor fields belong to cotangent bundle 
T*M, while their values are tensors related to base manifold M. If we replace T*M 
by tangent bundle TM, we can state another definition of extended tensor field. 

Definition 3.2. Extended tensor field X of type (r, s) in M is a tensor-valued 
function that maps each point q = (p, p) of some domain G C TM to a tensor of 
the space TJ(p, M). If G = TM, then X is called global extended tensor field. 

In the case of arbitrary smooth manifold M definitions 3.1 and 3.2 lead to differ¬ 
ent theories. But for Riemannian manifold M tangent bundle TM and cotangent 
bundle T*M are bound with each other by duality maps: 

g: TM T*M, g-i: T*M TM. (3.1) 

In local coordinates duality maps (3.1) are represented as index lowering and index 
raising procedures in arguments of extended tensor field X: 


n n 

Pi = ^ 9i] p’, p" = Y1 5 *^ Pj- 

i=i i=i 

Due to duality maps (3.1) two objects introduced by definitions 3.1 and 3.2 are 
the same in essential. We call them covariant and contravariant representations of 
extended tensor field X. 

In local coordinates x^, ... , x" extended tensor field X is represented by its 
components Xj); (a:\ ... ,a;”,pi,... ,p„) or Xj); *^(a:\ ... ,a:",p\... ,p"), de¬ 
pending on which representation (covariant or contravariant) is used. Extended 
tensor field X is called smooth if its components are smooth functions. 

Smooth extended tensor fields form a ring, we denote it by U = ^{T*M) in the 
case of covariant representation, and by j? = ^(TM) in the case of contravariant 
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representation. The whole set of smooth extended tensor fields in M is equipped 
with operations of 1) summation, 2) multiplications by scalars, 3) tensor product, 
4) contraction. It forms bi-graded algebra over the ring We denote this algebra 
by T(M) and call it an algebra of extended tensor fields in M: 


T(M) = 00^:w- (3-2) 

r=0 s=0 


Definition 3.3. A map D: T(M) ^ T(M) is called a differentiation of extended 
algebra of tensor fields, if the following conditions are fulfilled: 

(1) concordance with grading: D{Tf{M)) C Tf{M)] 

(2) R-linearity: D(X + Y) = L>(X) + D{Y) and D{XX) = XD(X) for A G R; 

(3) commutation with contractions: D{C(X.)) = C{D(X.)); 

(4) Leibniz rule: D(X (g) Y) = D(X) (g Y + X (g D(Y). 

Theory of differentiations in extended algebra of tensor fields (3.2) is considered in 
Chapters II, III, and IV of thesis [23]. In this section below we shall mention some 
facts from this theory needed for further use. 

Suppose that T(M) is extended algebra of tensor fields in M taken in contravari- 
ant representation. Then the set of its differentiations D{M) possesses the structure 
of module over the ring ^{TM). The set of extended vector fields (i. e. summand 
Tq{M) in direct sum (3.2)) also possesses the structure of ^{TM)-module. There¬ 
fore the following definition is consistent. 

Definition 2.1. Covariant differentiation V in the algebra of extended tensor 
fields T(M) is a homomorphism of ^{TM)-modules V : Tq{M) S)(M). Image 
of vector field Y under such homomorphism denoted by Vy is called covariant 
differentiation along vector field Y. 

For each covariant differentiation the expression VyX is lJ(rM)-linear with 
respect to Y. Therefore V can be treated as a map V : Tf{M) Tfj^ffM). 
Each smooth manifold M possesses exactly one canonical covariant differentiation 
V which is called vertical gradient. In local coordinates it is expressed by formula 






dp^ 


(3.3) 


In order to define other covariant differentiations one need some additional geo¬ 
metric structures in M. Thus, if M possesses affine connection T, one can define 
horizontal gradient V. In local coordinates it is expressed by formula 




a ^Zi... ir 


n n a 

a j^b 
qa 

a=l b—1 


dxi dp^ 

E \ ' pik _ \ ' \ ' pbfc vil . 

QO-k ji . js Z^ qjk jl-.-bk 


+ 


k—1 ak — 1 


k—1 bk — 1 


(3.4) 


If we take covariant representation of the algebra of extended tensor fields T(M), 
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then formulas for vertical and horizontal gradients are transformed as follows: 


_ ji--- js 

jl---js Qp^ 


(3.5) 




dx;i: 


dx^i 


a—1b—1 


dxi 


■Js 


dpb 


r*'“ 

Z_^ Qo^k ^^31. 
k—1 ak—1 


s n 

\ ^ \ ^ ybk vii . ir 

Qjk bk-.-Js ' 

k—1 bk — 1 


(3.6) 


In the case of arbitrary smooth manifold M gradients defined by formulas (3.3) 
and (3.4) are not related to those defined by formulas (3.5) and (3.6). However, if 
M is Riemannian manifold, then V and V defined by these two ways appear to be 
the same differentiations^ in different representations of algebra T(M). This fact 
is expressed by the following commutation relationships: 


V(Xog) = (VX)og, V(Xog) = (VX)og. 

Here g is duality map (3.1) defined by metric tensor of Riemannian manifold. 

4. Legendre transformation. 

Legendre transformation is usually used to relate Lagrangian and Hamiltonian 
dynamical systems. Suppose that M is smooth manifold and let be smooth 

extended scalar field in M taken in contravariant representation. Then dynamical 
system in tangent bundle TM described by differential equations 


dt \dvY dx^ 


is called Lagrangian dynamical system. Let’s apply covariant derivative (3.5) to L. 
As a result we get covector field p = p(p, v) with components 

dL 

= (4-2) 

If pair (p, v) is a point of tangent bundle TM , then pair (p, p) is a point of cotangent 
bundle T*M . This means that derivatives (4.2) determine a map 

X:TM -^T*M. (4.3) 


This map is known as Legendre transformation (see [30]). Below we assume Le¬ 
gendre transformation (4.3) to be invertible. Moreover we assume inverse map 

: T*M TM (4.4) 


to be smooth. Under these assumptions we can treat direct and inverse Legendre 
transformations (4.3) and (4.4) as nonlinear analogs of duality maps (3.1). Lagrange 


^One should only lower index q in (3.5). 
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function L{p,v) and Lagrange equations (4.1) are associated with tangent bundle 
TM. Vector v in arguments of Lagrange function is called velocity vector, while 
covector p with components (4.2) is called momentum covector. This gives rise to 
the following terminology. If X is an extended tensor field in contravariant repre¬ 
sentation and if Y = X»A“^, then we say that Y is p-representation or momentum 
representation for X, while X is called v-representation or velocity representation 
for Y. For the case of general smooth manifold M (without Riemannian metric) 
direct and inverse Legendre transformations (4.3) and (4.4) bind the following two 
representations of extended tensor fields: 


covariant p-representation 


contravariant v-representation 


If M is Riemannian manifold, we have four representations per each extended field: 


covariant p-representation 

contravariant v-representation 

contravariant p-representation 

covariant v-representation 


Now let’s consider the following two extended scalar fields h and H: 

n 

h = '^v^ViL-L, H = hoX-\ (4.5) 

i=l 

Scalar field H is known as Hamilton function, while h is its v-representation. La¬ 
grange function L and its p-representation I = L o can be expressed through 
Hamilton function H by formulas similar to (4.5): 


l = ^p,V^H - H, L = loX. (4.6) 

i=l 

Applying Legendre transformation to Lagrangian dynamical system (4.1), we can 
transform them to Hamiltonian dynamical system in cotangent bundle T*M: 

, dH dH 

^ 

opi ox^ 

This fact is well known (see [30]), as well as above formulas (4.5) and (4.6). 

5. Modified Hamiltonian and Lagrangian 
DYNAMIGAL SYSTEMS IN RIEMANNIAN MANIFOLDS. 

Now suppose that manifold M is equipped with Riemannian metric g. All results 
we discussed in section 2 were obtained under this assumption. Remember that 
wave front dynamics is described by modified Hamilton equations (2.8). Using 
(3.5) and (3.6), we can replace partial derivatives in them by covariant derivatives: 




VtPi 




X' 


n 


n 


(5.1) 
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Time derivatives , ... , i" in (5.1) are components of tangent vector to trajectory 
p = p{t), while Vt is standard covariant derivative with respect to parameter t 
along this curve. Similar to original Hamilton equations (4.5), modified Hamilton 
equations are associated with cotangent bundle T*M. Using inverse Legendre map 
(4.4), one can transform them to v-representation. This was done in paper [4]. As 
a result modified Lagrange equations were obtained: 


V.(V.L) = 

Denominator D in original p-representation is given by formula 

n 

Upon passing to v-representation in (5.2) formula (5.3) transforms to 


(5.2) 


(5.3) 


0 = (5.4) 

i=l 

Formula (5.4) means that we can express D in v-representation explicitly through 
Lagrange function L. 


6. Newtonian dynamical systems admitting 

NORMAL SHIFT IN RiEMANNIAN MANIFOLDS. 

One can see that vector v for modified Lagrangian dynamics (5.2) do not coincide 
with actual velocity vector. If we denote actual velocity vector by u, then we derive 


u 


i 


n’ 


Formula (6.1) defines nonlinear map similar to Legendre map A: 


( 6 . 1 ) 


p-.TM^TM (6.2) 

If this map (6.2) is invertible and if inverse map is smooth, then one can 
transform modified Lagrange equations (5.2) to the following form: 


i* = u*, Vtu* = U*(x^,... , a;", ... , m"). (6.3) 

Differential equations (6.3) determine Newtonian dynamical system in M. Ex¬ 
tended vector field F with components F^, ... , F” in (6.3) is called force field of 
this Newtonian dynamical system. 

Note that nonlinear map (6.2) is more complicated than Legendre map. Almost 
each modified Lagrangian dynamical system can be transformed to Newtonian form 
(at least locally). However, converse is not true. Moreover, even if it is known that 
Newtonian dynamical system (6.3) is derived from modified Lagrangian dynamical 
system (5.2), there is no explicit formula for F. In paper [4] one can find explicit 
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formula for F, but in very special case, when Lagrange function L is fiberwise 
spherically symmetric with respect to Riemannian metric g: 

Fk = -H-^^^-{2N^Nk-5i). (6.4) 

Here W = W{p,\vl\) = h° is u-representation for Hamilton function H and W 
is partial derivative of W with respect to its (n + l)-th argument u = |u|, fV* and 
Nk are components of unit vector N = u/|u|. 

Remember that in the example considered in section 2 Hamilton function H is 
given by formula (2.9). It is fiberwise spherically symmetric, i.e.H = H(p, |p|). 
Applying inverse Legendre map we get fiberwise spherically symmetric function 
h = h{p, |v|) = H o X. Its u representation then is fiberwise spherically symmetric 
function W in formula (6.4). Thus we have the following theorem. 

Theorem 6.1. For generalized wave equation (2.3) in Riemannian manifold M 
wave front dynamics in the limit of short waves is described by Newtonian dynamical 
system (6.3) with force field (6.4). 

Theorem 6.1 is the main result of paper [4]. It establishes a link between wave 
propagation phenomena and the theory of dynamical systems admitting normal shift 
(see papers [7-22]). Below we give brief introduction to this theory. 

Let M be Riemannian manifold and let a be some smooth hypersurface in M . 
Suppose that p is a point of a and n(p) is a unit normal vector to a at the point p. 
Under these assumptions we can consider initial data 

xM =x^{p), mM =v{p)-n'^{p) (6.5) 

li=0 lt=o 

for Newtonian dynamical system (6.3). Similar to initial data (2.10) for modihed 
Hamiltonian dynamical system, here initial data (6.5) define a shift of hypersurface 
a along trajectories of Newtonian dynamical system (6.3). This shift is called nor¬ 
mal shift if hypersurfaces ctj, which are obtained from a by shift, keep orthogonality 
to shift trajectories in time. 

Definition 6.1. Newtonian dynamical system (6.3) is called a system admitting 
normal shift if for any hypersurface cr there is a smooth function v = v{p) on a such 
that initial data (6.5) with this function v define normal shift of cr along trajectories 
of dynamical system (6.3). 

Suppose that po is some fixed point of hypersurface a and let vq be some fixed 
constant. Let’s normalize v{p) by the following condition: 


J^(po) = ^ 0 - (6.6) 

Definition 6.2. Say that Newtonian dynamical system (6.3) satisfies strong nor¬ 
mality condition if for any hypersurface cr, for any point po G cr, and for any constant 
t'o 7 ^ 0 there is a smooth function v = i'{p) on a normalized by the condition (6.6) 
and such that initial data (6.5) with this function v define normal shift of cr along 
trajectories of dynamical system (6.3). 

Strong normality condition, in contrast to the normality condition from defini¬ 
tion 6.1, is less obvious. But it is more convenient for to study by mathematical 



14 


RUSLAN A. SHARIPOV 


methods. In papers [12] and [13] the following two systems of differential equations 
for the force field F of Newtonian dynamical system (6.3) were derived: 


E + Pfc=0, 


i=i 


E E (^*^1 + - 2 Fj) PI + 

i=i j=i 

+ E E (^ - E f\ pi = 0, 

i=l j—\ \ r—1 / 


(6.7) 


^ ^ ^ Z? Y7 fT 

EE^-E - v,F, 1 = 


i=i j=i 


V m—l 




i=i j=i 


n n n 


E E = E E E 


L i=ii=i 


j-'=l m=l 


n — 1 


( 6 . 8 ) 


The equations (6.7) were called weak normality equations, while other equations 
(6.8) were called additional normality equations. In Chapter V of thesis [23] the 
following theorem was proved. 

Theorem 6.2. Newtonian dynamical system (6.3) satisfies strong normality condi¬ 
tion if and only if its force field F satisfies complete system of normality equations 
consisting of weak normality equations (6.7) and including additional normality 
equations (6.8) in the case of higher dimensions n ^ 3. 

Weak normality equations (6.7) are related to weak normality condition. In order 
to formulate this condition let’s consider one-parametric family of trajectories of 
Newtonian dynamical system (6.3). Denote it as follows: 


p = pit,y)- 


(6.9) 


Here t is time variable and j/ is a parameter. In local coordinates this one-parametric 
family of trajectories (6.9) is expressed by functions 


= x'^{t,y), 


x” = x^{t, y). 


( 6 . 10 ) 


Differentiating (6.10) with respect to parameter y, we get vector t with components 
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Vector T is called vector of variation of trajectories. Then from (6.3) we derive 


n n n 

9=1 i=l i=l 
n n 

+ VtT^ VqF^ + 

9=1 9=1 


( 6 . 12 ) 


Here components of curvature tensor for metric g, while and are com¬ 
ponents of velocity vector u. Function ip defined as scalar product 

n n 

^p = {u\T) = Y'^9^Ju"r^ (6.13) 

i=i i=i 

is called function of deviation. From (6.12) one can derive the following ordinary 
differential equation for the function of deviation (6.13): 


2n 

^Q(t)(p«=0. (6.14) 

i=0 

For general Newtonian dynamical system (6.3) this is homogeneous ordinary differ¬ 
ential equation of the order 2n (see details in Chapter V of thesis [23]). However, 
in special cases the equation (6.14) can reduce to lower order differential equation. 
Weak normality condition below specifies one of such cases. Indeed, let’s consider 
some trajectory p = p(t) of Newtonian dynamical system (6.3). It can be included 
into one-parametric family of trajectories (6.9) by various ways. This defines vari¬ 
ous variation vectors r with components satisfying differential equations (6.12) and 
various deviation functions (6.13) on the trajectory p = p{t). 

Definition 6.3. Say that Newtonian dynamical system (6.3) satisfies weak nor¬ 
mality condition if for each its trajectory p = p{t) and for any vector of variation t 
on this trajectory corresponding function of deviation ip(t) satisfies homogeneous 
second order ordinary differential equation 


(p = A{t) (fiB{t) ip (6.15) 

with coefficients depending only on choice of trajectory p = p{t). 

As it was shown in paper [12], weak normality condition is equivalent to weak 
normality equations (6.7) for the force field of Newtonian dynamical system (6.3). 

Now let’s proceed with additional normality condition. In order to formulate this 
condition let’s consider some smooth hypersurface a in M and let’s fix some point 
Po on a. Denote by ... , local coordinates on a in some neighborhood of 
fixed point po- Setting up initial data (6.5), we can define a family of trajectories 
p = p{t, y^,. ■. , p") of Newtonian dynamical system (6.3) starting at the points of 
a. Now this is (n — l)-parametric family of trajectories expressed by functions 

( =a;i(t,p\... ,p""^). 


,p” ^)- 


(6.16) 
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in local coordinates x^, ... , x” in M. Differentiating these functions (6.16) with 
respect to parameters y^, ... , as it is done in (6.11), we get n — 1 variation 

vectors Ti, ... , t„_i. It’s easy to note that variation vectors Ti, ... , t„_i form a 
frame of tangent vectors to hypersurfaces ctj obtained by shifting initial hypersurface 
a along trajectories (6.16) of dynamical system (6.3). Therefore corresponding 
deviation functions (/?i, ... , ipn-i serve as measure of orthogonality of Ut and shift 
trajectories (6.16). They should be identically zero in order to provide orthogonality 
of shift: .. , y"“^) = 0. If we consider initial data 


ipk\ =0, i^fcl =0, (6.17) 

li=0 lt=0 

then we can see that first part of initial conditions (16.17) is fulfilled due to ini¬ 
tial data (6.5). Second part of these conditions can be transformed to differential 
equations for the function v = v(p) = i'{y ^,... , y”“^) in (6.5): 

F)ij 

^ = -^-i(F|x,). (6.18) 

If dimM = n ^ 3, then the equations (6.18) form complete system of Pfaff equa¬ 
tions for scalar function v. The condition of its compatibility is known as additional 
normality condition. 

Definition 6.4. Say that Newtonian dynamical system (6.3) satisfies additional 
normality condition if for any smooth hypersurface a in M and for any local coor¬ 
dinates 2 /^, ... , on cr corresponding Pfaff equations (6.18) are compatible. 

In paper [13] it was shown that for n ^ 3 additional normality condition is 
equivalent to additional normality equations (6.8) for the force field of Newtonian 
dynamical system (6.3). In two-dimensional case n = 2 situation is quite different. 
Here we have only one parameter y = y^ and (6.18) turns to unique ordinary 
differential equation, which is compatible with itself in anyway. Therefore in two- 
dimensional case additional normality condition is always fulfilled. This special 
case is studied in thesis [24]. 

In higher dimensional case n ^ 3 complete system of normality equations in¬ 
cludes both (6.7) and (6.8). For this case in Chapter VII of thesis [23] explicit 
formula for general solution of complete system of normality equations was derived: 


Fk 


hjW) Nk 
W' 


H-i: 

i=l 


V,IT 

W' 


(2N^Nk 


SI)- 


(6.19) 


Comparing formulas (6.4) and (6.19) we see that they are quite similar. They differ 
only by first term in (6.19), where h = h{w) is arbitrary function of one variable. 
This fact indicates that modified Lagrangian dynamical systems (5.2) describing 
wave front dynamics and Newtonian dynamical systems (6.3) admitting normal 
shift of hypersurfaces in Riemannian manifolds are closely related with each other. 
In further sections we are going to reveal this relation in more general case, when 
manifold M is not equipped with Riemannian metric. Problem of interpreting first 
term in (6.19) should be considered in separate paper. 
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7. Weak normality phenomenon for 
MODIFIED LAGRANGIAN DYNAMICAL SYSTEMS. 

Let M be a smooth manifold which is not equipped with Riemannian metric, but 
which is equipped with modified Lagrangian dynamical system. Let L = L(p, v) be 
Lagrange function for this system. This is extended scalar field in contravariant v- 
representation. In the absence of Riemannian metric we cannot use spatial gradient 
(3.4). Therefore we write modified Lagrange equations as 


X 


i 




1 dL 

dt \ dv^) fl dx^' 


Denominator D in (7.1) is determined by formula (5.4): 






i=l 


dL 

dv^ 


(7.1) 


(7.2) 


Remember that formula (4.2) defines Legendre map (4.3). Below we assume this 
map A to be invertible, and moreover, we assume inverse map to be smooth. 
Local invertibility of A means that matrix g with components 


Vi VtL 
Lij = —^ 


1 d^L 

2 dv'^ dv^ 


is non-degenerate. We shall assume this matrix to be positive: 


(7.3) 


/X > 0. 


(7.4) 


For real mechanical systems this condition /x > 0 is fulfilled since kinetic energy K 
of such systems is positive quadratic function of velocity vector. For such systems 
components of matrix /x do not depend on v, hence one can choose /x to be metric 
tensor for Riemannian metric in M. However, we shall consider more general case, 
when fjb is extended tensor field with components depending on v. 

In addition to inequality (7.4) we shall assume that denominator Q in modified 
Lagrange equations (7.1) (which is determined by (7.2)) is positive function: 

D > 0, for V 7 ^ 0. (7.5) 


This assumption is consistent since for real mechanical systems Lt = 2K. 

Now let’s consider one-parametric family of trajectories p = p{t, y) of modified 
Lagrangian dynamical system (7.1). In local coordinates these curves are expressed 
by functions (6.10). Formula (6.11) then defines vector of variation t. In order 
to define function of deviation p we could use formula (6.13) with matrix (7.3) as 
metric. However, we choose another formula for p: 

n 

<L= {p\t-) = C^-S) 

i=l 

Here p is momentum covector defined by formula (4.2), while angular brackets 
denote contraction of vector and covector^. Function (p in (7.1) can be treated as 


^Such notations are often used in quantum mechanics. See [31]. 
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scalar product of vectors v and t. This scalar product is linear with respect to 
vector T, but it is nonlinear respect to vector v. Such scalar products usually arise 
in Finslerian geometry (see Chapter VIII of thesis [23]). 

Let V = v(t, y) be velocity vector for one-parametric family of trajectories p = 
p{t,y) of modified Lagrangian dynamical system (7.1). In local coordinates this 
vector-function is expressed by the following scalar functions: 


v'^ = v'^{t,y), 




(7.7) 


Differentiating (7.7) with respect to parameter y, we get series of functions 



(7.8) 


In contrast to r* in (6.II), these functions 6*^, ... , 0" are not interpreted as com¬ 
ponents of vector. We shall use them in order to simplify further calculations. 
Differentiating first equation (7.1) with respect to y, we obtain 


Q ^ dv^ 


n n 


-EE 


d^L 

dyk QyS 


yi yk QS 

D2 


n n 


EE 


d^L 

dyk Qj.S 


yt yk yS 


(7.9) 


Differentiating second equation (7.1) with respect to parameter y, we get 


E 


dyt QyS 


n „2r « " 

+EE 


d^L 


n n 


k—1s—1 
n n 


Qyl QyS Qyk 


r + 


o3r " " o3r 

V V_ — _-L V y_ — _r"' -b 

f),,i f)„s f)rrk „ r, U T -r 


k—1 s—1 


k—1 s—1 


dv'- dx^ dv^ 


EE 

k—1 s—1 


d^L 


dv^ dx‘^ dx^ 




_ " ^ Q^L 


n dx^ dv^ 


^ I d^L 

E o " ^ 


S = 1 


n dx^ dx^ 


^ dL dL 

2-^ rin^s Q2 


s=l 
n n 


dL d% 9’^ V y 

Z^ Z^ Q^i Qyk QyS Z^ Z^ Q^i Qyk Q^s 

k—1 s=l k—1 s—1 


(7.10) 


Both (7.9) and (7.10) form a system of homogeneous linear ordinary differential 
equations with respect to functions r^, ... , t” and 9^, ... , 0”. This system of 
equations is an analog of equations (6.12) considered above. 
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Function of deviation Lp defined by formula (7.6) depends linearly on components 
of vector t. Let’s calculate time derivatives of this function. For ip we get 


. ^ dL ^ d ( dL 

S = 1 ^ 




9a;'* n 9w® fl dv^ dv^ 

s—1 s=l s—1 r—1 


n n n 


sr^sr^sr^ dL d‘^L v>^ 9^ y- y" y- d‘^L 

A-^ QyS Qyk Qyr fl 2 A-^ 2 -^ A-^ QyS Qyk Q^r Q 2 

S—1 k—1 r=l s—1 k—1 r—1 


dL T« ^ 9L r'* d^L v’^ 9^ d^L 


E UJ-j T V UlJ T V V 

r)nrS O r ^ O ^ ^ ^ 


dx^ n 9a;'* O 9u* dv^ id ^ ^ dv^ 9a;’’ il 

S—1 s=l k—1 r=l k—1 r=l 


-EE- 


V T 


In the above calculations we used second equation (7.1), used formula (7.2) for 
denominator il, and used differential equations (7.9). As a result we obtained 


9L T® d'^L A A d% 


. V—^ iJJ-J I ^^> 

^ O ' 


dx^ Q ^ ^ dv^ dv^ ^ ^ dv^ dx^ Q 

s—1 k—1 r—1 k—1 r=l 


-EE 


V T 


(7.11) 


Now let’s differentiate (7.11) once more. This yields 


" " o2r ns n n 

p = -yy ^-yy 

f)^}^ O f ^ f ^ 


dv'^ dv^ ^ ^ ^ dv'^ dx^ ^7 ^ dx^ ^ 

i—l s—1 i—1 s—1 s—1 


’V" —— 

f ^ f)nr^ O 


n n n 


^ ^ ^ d^L v'^v^e^ 

2 ^ 2 -^ 2 -^ Qyk Qyl QyS ^ 2 -^ 2 -^ 2 -^ Qj,k Qyt QyS ^ 

k—1 i—l s—1 k—1 i—l s—1 


n n n 

-EEE 

k—1 i—l s=l 


r^'^-r • Zc 7 c 77 n n 

QiL yk yt 


Qyk Qyl Q^S ^ 


-EEE 


d^L x'^ w’ T® 


k—1i—l s—1 


dx^ dv^ 9a;'* il 


2-^ 2-^ f)q)i f)q)S Q 2-^ 2-^ 


d'^L fi* r'* 


Qyl QyS Q A-^ A-^ Qyl Qj.S ^ 

i—l S=1 7=1 s—1 


-EE- 


d^L V 


V T 


EE- 


d^L i'*’ il 


dv^ dx^ fl ^ ^ dx^ dx^ il 

k—1 s—1 k—1 s—1 


Using equations (7.10), we can eliminate all entries of derivatives 0® from the above 
expression for i)3. As a result we get reduced formula 


7^ = -EE 


d^L 

9a;’ 9u® 



EE 

i=l s=l 


d^L u’ T® 
9a;’ 9x'* 
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dL dL 

2-^ 2-^ rtqr^ r)i^s 


n n n 


V W — 

' pirni 


9 x* i 9 w® ^ ^ ^ dx^ dv^ dv^ Vfi 

i—1 s—1 i—1 k—1 s=l 


n n n 


W V ‘^L d^L v^ v^ r" ^ dL f" 

^ ^ ^ /9t** F)n^ Pi'T^ O^ '' ^ 


z=l fc=ls=l 


Sx'* Q 


(9t>* 9?;® ri ^ (9a;^ 9x® O 

z=l S—1 k—1 s=l 

Now we eliminate all entries of f®, using equations (7.9) for this purpose: 


d^L v ^ T® 


Tl Tl T ? /I Q TT' TT' 

9x* 9n® i9a;* dx^ 

2=1 S=1 2=1 S = 1 


Vt 

9x® ^ 9n* 9w® n ' ^ ^ dx^ 9a;® H n 

S=1 2=1 S=1 fc = l s=l 


uL U ^-~\ ^-~\ 

2-^ pr^S 02 ~ 2-^ 2-^ 


EE 


‘P- 


(7.12) 


Then, using first equation in (7.1), we express time derivative x^ through . As a 
result formula (7.12) for (p reduces to the following one: 


.. _ d^L n* 0® ^ dL 6»® 

^ 9x* 9n® ^ 9x® 

2=1 S=1 S=1 


-EE 

2=1 S—1 


d‘^L V* 6»® 


(7.13) 


9n* 9w® n n 


7^- 


Now, if one take into account second equation (7.1) written in expanded form, then 
three terms in right hand side of (7.13) can be canceled. This yields 


(p+^ip = {). (7.14) 

Thus, in the end of huge calculations we get very simple relationship (7.14), which 
is homogeneous second order linear ordinary differential equation. It is even more 
simple than analogous equation (6.15) considered in previous section. Now, if we 
formulate definition 6.3 respective to modified Lagrangian dynamical system (7.1) 
and if we use formula (7.6) for (/?, then from (7.14) we derive the following theorem. 

Theorem 7.1. Each modified Lagrangian dynamical system (7.1) satisfies weak 
normality condition with respect to deviation functions (7.6) determined by its 
own Lagrange function L. 


8. Additional normality phenomenon. 

In order to reproduce results of section 6 in present more complicated geometric 
environment we should consider some hypersurface cr in M, and we should arrange 
a shift of cr by means of modified Lagrangian dynamical system (7.1). Fortunately 
we should not invent something absolutely new for this purpose. Wave front dy¬ 
namics considered in section 2 suggests a way of how to do this. In the absence of 
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Riemannian metric we cannot choose unit normal vector on a. However, we can 
take normal covector n = n(p), which is unique up to a scalar factor. Then we can 
set up Cauchy problem with the following initial data: 


xM =x\p), pi\ =v{p)-ni{p) (8.1) 

li=0 lt=0 

(compare with (2.10) above). Here p is a point of cr and Pi are components of 
momentum covector p defined by formula (4.2): 

^ . dL 

”• = = a?- 

Initial data (8.1) determine initial velocity v implicitly through initial momentum 
covector p due to invertibility of Legendre map A. Applying initial data (8.1) to 
modified Lagrangian dynamical system (7.1), we obtain a family of trajectories 
p = p(t,y^,... ,y") starting at the points of hypersurface a. Similar to (6.16), in 
local coordinates these trajectories are expressed by the following functions: 


=x^(t,y\... ,y” 1), 


= x^{t,y\... ,y” 1). 


( 8 . 2 ) 


These functions (8.2) determine variation vectors ti, ... , t„_i with components 

i _ 

1 Qyj 

(compare with (6.11) and (7.8)). Each variation vector determines correspond¬ 
ing deviation function according to the formula (7.6). We denote these deviation 
functions by pi, ... , p„_i as in section 6. 

Definition 8.1. Shift of initial hypersurface a determined by modified Lagrangian 
dynamical system (7.1) and by initial data (8.1) for it is called normal shift in inner 
geometry of dynamical system (7.1) if all deviation functions pi, ... , pn-i are 
identically zero. 

Due to differential equation (7.14) for deviation functions in order to arrange a 
normal shift of cr it is sufficient to provide initial conditions 


Pi, I =0, pd = 0 (8.3) 

lt=o lt=o 

just the same as in (6.17). First part of initial conditions (8.3) is fulHlled due 
to initial data (8.1). Second part of these conditions should be transformed to 
differential equations for the function n = ^(p) = v{y^, ■ ■ ■ , y”~^) in (8.1). For this 
purpose we could use formula (7.11) derived in section 7. However, initial data 
(8.1) explicitly relate function n = v{p) with initial value of momentum covector 
p, while relation to velocity vector v is implicit. Therefore it is easier to transform 
formula (7.6) to p-representation. This yields 

n 

‘Pr = ( p | t *) = ^PsTf. 


(8.4) 
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Remember that modified Lagrange equations, when transformed to p-representa- 
tion, look like modified Hamilton equations (2.8): 


1 dH 
H dpi ’ 


1 dH 
H dx'‘ 


(8.5) 


Here Hamilton function H is determined by formula (4.5) and H is given by formula 
(5.3). Now, differentiating formula (8.4), we obtain 


n n n 

Vi=^PsTt = ^PsT^ + '^Ps t / 

S=1 S=1 S=1 


E 


1 dH ^ 




( 8 . 6 ) 


In order to calculate time derivatives ff in formula (8.6) we use first part of modihed 
Hamilton equations (8.5). A a result for ff we get 



Let’s substitute this expression into (8.6). This yields 


dH ^ ^ ^ ^ d f I dH 

= &7(n^ 

s=l r—1 ^ 


EE 


Ps 


dpr d f 1 dH 


+E 


5=1 r—1 

dpr d 


dy^ dpr <9;?, 


= E' 

r=l 


d 

dx'^ 


^ fl dps 

<S—1 ^ / 


^Ps dH\ dH dps A 1 dH 

O drj. 1 .O dr), d?/* ^ , 


^ dy^ dpr^j^^^ dps j ^ dps dy^ 


S = 1 


H dx'^ 


First two terms in right hand side of the above equality are identically zero. This 
follows from formula (5.3) for Q. Thus for (pi we have 


^ H dps dy^ ^ H dx'^ 

If we recall initial conditions (8.3), then from (8.7) we derive 


(8.7) 


s + 

^ dx^ ^ dps \dy'^ 


S = 1 


= 0 . 


t=o 


( 8 . 8 ) 


Calculating partial derivatives dps/dy'^ in (8.8), we should remember (8.1). Then 

(8.9) 


dps 

dy’’ 


dv dus 1 du dus 

t=o dy^ dy^ v dy'' ^ dy"^ 
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Substituting this expression into (8.8) and using formula (5.3) for we can trans¬ 
form (8.8) to the partial differential equations for v. 

= ——-V —^ 18 101 

jy Qyl Qyl Qp ^ Qy-S 

Differential equations (8.10) are analogs of the equations (6.18). If n ^ 3, then 
they form complete system of Pfaff equations for the function v = v(y ^,... , ?/"“^). 
Therefore we can formulate additional normality condition for modified Lagrangian 
dynamical system (7.1) as compatibility condition for Pfaff equations (8.10). 

Suppose that n ^ 3. Let’s examine if differential equations (8.10) are compatible. 
For this purpose let’s calculate second order partial derivatives of v using (8.10): 


d'^v 
dy^ dyJ 


^ dH d'^Us 
“ D dps dy’- dy^ 


E 


dH V 

dx^ Q dy^ dyi 


Vfi dps dpr dy^ dy^ dx^ dx^ * ^ 

s—1 r—1 ^ ^ g—i r—1 


^^ dH dH dus 

2^ 2^ g^s gp gyj * 2^ 2^ gp g^r gyi 

S—1 r—1 s—1 r—1 


- 


^ ^ 2 ^ f ^ \ ^Ps dUr ^ ^ 9/1 dH \ dps r 

dPs \ ^ ^Pr ) dy’- dyi ^ dps V D dx^ J dy^ 


*-VV —(22IL\ ^ r 

“ “ dx’^ \ D dpr J dy^ ^ ^ dx^ V D dx"^ J 


For Pfaff equations (8.10) to be compatible, right hand side of the above equal¬ 
ity should be symmetric in indices i and j. First four terms there are obviously 
symmetric. Below we shall not write such terms explicitly denoting them by dots: 


d'^v 
dy^ dy^ 


n n / 

+ee( 


2v^ dH dH 
dx‘^ dpr 


d 

dx^ 




_d_[ 2 2E\ 2L2h: 

^Ps dpr) dy^ dyi dx^J dy^ 

E E gp I fi gp J gyt gyj E E g^s I ^ g^r J 

s=l r—1 ^ L / >j J r—1 ^ 


In the above calculations we used formula (8.9) for partial derivatives dpsjdy^. 
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Below we use the equations (8.10) for to express partial derivatives dvjdif-. 
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dy'^ dyi 


n n 

EE 
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d fl dH 


s=l r—1 

2 - d fl dH\ dH\ dn. 


+ 


E fl Qp^ yfi dprj dx^ J dy^ ’* 


-rf 


+EE 


dx^ \ Q dp, 

' dH dH 
dps dx'^ 


d fl dH 


Pa d /I dH \ dH \ dn, 


2 u I L un \ V ly Pq u ! L un \ un uiLg r 

fin. I O /9'r^ / f ^ O fin. , O fi'r'^ ) fin. fiif ^ 


dps V ^ dx^ J “ Q dpq \ dx'^ J dps j dy 


EE E 

s—l r=l \ q—1 


Pq d f 1 dH \ dH jy^ dQ dH \ dris dn^ 


dpq \ dpr J dps dps dpr I dy^ dy^ 


n dpq \ n dx^ J dx^ t^2 fl^rS rtrrr 


S—l r—1 \ q—1 


02 dx^ dx^ I ■ 


Now we are able to write compatibility condition for Pfaff equations (8.10). It 
breaks into three separate parts. These are the following equalities: 


21/2 dH dH 
n2 i9a;® dpr 

_ 1/2 dH dH 

~'^dfdf 


d n dH 


E 


Pq d fl dH\ dH 


dx^ \^l dprJ dpq \ Q dprJ dx 


d fl dH 


dpr dx^ 


E' 

9=1 


Pq d fl dH\ dH 


fl dpq \ fl (9x® J dpr ’ 


^ v^Pq d fl dH\ dH 1/3 dn dH _ 

“ n dpq \n dpr) dps 112 Qp^ Qp^ 


= E 


ly^Pq d fl dH\ dH ^ dfl dH 


^ fl dpq \ fl dpsj dpr 112 Qp^ Qp^ 


^vpq dfl dH\ dH V dfl dH _ 

H dpq \ H dx'' J dx'' 112 Qj. 8 Q^r 


_^vpq d fl dH\ dH ly dfl dH 

“ n dpq \11 dx" J dx'" H2 dx'' dx" 


( 8 . 11 ) 


( 8 . 12 ) 


(8.13) 


It is easy to check that all these three equalities (8.11), (8.12), and (8.13) turn to 
identities if we substitute (5.3) for H. Therefore we can formulate the following 
main result of this section. 

Theorem 8.1. Each modified Lagrangian dynamical system (7.1) satisfies addi¬ 
tional normality condition with respect to deviation functions (7.6) determined 
by its own Lagrange function L. 
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One should note here that PfafF equations (8.10) are not only compatible, but 
they are also explicitly integrable in form of equality H{p,v ■ n(p)) = const, which 
is similar to the equality (2.11). 

9. Summary and conclusions. 

Theorems 7.1 and 8.1 form main result of present paper. Now we are to under¬ 
stand this result. Thus, we have arbitrary smooth manifold M without Riemannian 
metric, but equipped with Lagrange function L = L{p,'v) defining invertible Le¬ 
gendre transformation A and satisfying two conditions (7.4) and (7.5). It’s clear that 
these conditions are rather non-restrictive. Despite to the absence of Riemannian 
metric, under the above assumptions 

1) one can define concept of normal shift with respect to geometric structures 
determined only by Lagrange function L; 

2) one can formulate weak and additional normality conditions] 

3) one can prove that modified Lagrangian dynamical system (7.1) satisfies both 
normality conditions with respect to geometric structures determined by its 
own Lagrange function L; 

The results listed above generalize a part of theory of dynamical systems admitting 
normal shift from Riemannian geometry to the geometry of Lagrangian dynamics. 
However, this is not complete generalization. Indeed, in Riemannian case geometry 
was determined by metric tensor g, while dynamics was determined by force field F. 
Theory is based on the interplay of these two structures. Here both geometry and 
dynamics are determined by Lagrange function L yet. In further generalizations 
one should introduce another dynamical system in M (either Lagrangian or not 
Lagrangian), and then one should measure its capability to implement normal shift 
of hypersurfaces in geometry determined by L. 
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